Effects of ground-state correlations on collective excitations of °0 
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The effects of the correlations in the ground state of ^''O on the octupole and dipole excitations are 
studied using the extended random phase approximation (ERPA) derived from the time-dependent 
density-matrix theory. It is found that the ground-state correlation effects are significant especially 
in the octupole excitation. It is shown that the first 3~ state calculated in the random phase 
approximation (RPA) is shifted upward when the self-energy contributions are included in particle 
- hole pairs. The coupling to the two particle - two hole states plays a role in shifting the first 
3~ state down to the right position. It is also found that the dipole strength is fragmented due to 
the partial occupation of the single-particle states and that the peak position of the giant dipole 
resonance calculated in ERPA is little changed from that in RPA due to the above-mentioned 
competing effects: the increase in particle - hole energy and the coupling to two particle - two hole 
configurations. 

PACS numbers: 21.60.Jz, 21.10.Pc, 27.20.-f-n 



The random phase approximation (RPA) based on the 
Hartree-Fock (HF) ground state has extensively been 
used to study nuclear collective excitations. It is gener- 
ally considered that the HF -I- RPA approach is the most 
appropriate for double-closed shell nuclei for which the 
HF theory would give a good description of the ground 
states. However, recent theoretical studies for ^^O indi- 
cate that the ground state of ^^O is a highly correlated 
state [l|, 01 , which indicates the necessity of using beyond 
RPA theories to study collective excitations of ^^O. In 
this paper we investigate how octupole and dipole exci- 
tations of ^^O are affected by ground-state correlations, 
using an extended RPA (ERPA) that has been derived 
from the time-dependent density-matrix theory (TDDM) 
d, 13 • We show that the octupole excitation of ^^O is 
quite sensitive to the ground-state correlation effects. 

The TDDM consists of the coupled equations of motion 
for the one-body density matrix naa' (the occupation 
matrix) and the correlated part of the two-body density 
matrix Capa'P' (the correlation matrix). These matrices 
are defined as 

«aa'(t) = (*(0|aLaa|*(i)), (1) 

Ca/ia'fi'it) = {<i>{t)\al,al,afiac\<^{t)) 

- inaa'{t)nfjp'{t) - nap'{t)nap'{t)), (2) 

where |<I'(i)) is the time-dependent total wavefunction 
$(i)) = exp[-iHt]\^{t = 0)). The equations of mo- 
tion for reduced density matrices form a chain of cou- 
pled equations known as the Bogoliubov-Born-Green- 
Kirkwood-Yvon (BBGKY) hierarchy. In TDDM the 
BBGKY hierarchy is truncated by replacing a three-body 
density matrix with anti-symmetrized products of the 
one-body and two-body density matrices. The TDDM 
equation for Ca/ja'P' contains all effects of two-body cor- 
relations; particle - particle, hole - hole and particle - hole 
correlations. The ground state in TDDM is given as a 
stationary solution of the TDDM equations. The station- 
ary solution can be obtained using the gradient method 
[ij. This method is also used in the present work. The 



ERPA equations used here are derived as the small am- 
plitude limit of TDDM and are written in matrix form for 
the one-body and two-body amplitudes x'^^, and X^g^,g, 
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The one-body sector of Eq. ^ Ax^ — w^S'ix^ is for- 
mally the same as the equation in the self-consistent 
RPA (SCRPA) [B!], which includes the effects of ground- 
state correlations through Uaa' and Ca^a'p'- We refer to 
the approximation Ax^^ = tu^Six'^ as the modified RPA 
(mRPA). To explain the role of the correlation matrix in 
the mRPA equation, we explicitly show the matrices A 
and 5*1! 

A ( aa' : AA') = (e^ - eQ')("-a'a' ~ «qq)<^qA<5q'A' 
+ («Q'Q' - naa){nyx' - nxx){a\' \v\a' X) A 
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where the subscript A means that the corresponding ma- 
trix is antisymmetrized and Uaa' is assumed to be di- 
agonal. The first two terms on the right-hand side of 
Eq. ^ are the same as those in the RPA equation, the 
next two terms with Ca^a'fi' describe the self-energy of 
the particle - hole (p-h) state due to ground-state cor- 
relations [5], and the last four terms with Capa'p' niay 
be interpreted as the modification of the p-h interaction 
caused by ground-state correlations Q. The difference 



between mRPA and SCRPA is in the method to deter- 
mine n^a' and Ca^a'p''- in SCRPA these matrices are 
self-consistently determined using the amphtudes x'^^, 
whereas here we calculate them in TDDM. If the corre- 
lated ground state is replaced by the HF ground state, 
mRPA and ERPA become the same as RPA and the sec- 
ond RPA (SRPA) [6j], respectively. 

The occupation probability Uaa and the correlation 
matrix Capa'/S' are calculated within TDDM using the 
lP3/2i lPi/2i l'^5/2 and 2si/2 states for both protons and 
neutrons. For the calculations of the single-particle states 
we use the Skyrme III force. To reduce the dimension 
size, we only consider the two particle - two hole (2p-2h) 
and 2h-2p elements of Capa'P'- A simplified interaction 
which contains only the to and t^ terms of the Skyrme III 
force is used as the residual interaction. The spin-orbit 
force and Coulomb interaction are also omitted from the 
residual interaction. To avoid a cumbersome treatment of 
the rearrangement effects of a density-dependent force in 
extended RPA theories 7] , we use the three-body version 
of the Skyrme interaction, v^ = tj,5^{ri — r2)8'^{ri — rs), 
which gives the following density-dependent two-body 
residual interaction: tzPn5^{r — r'), t3p6^{r — r')/2 and 
t3PpS^{r — r') for the proton-proton, proton- neutron and 
neutron- neutron interactions, respectively, where pp, pn 
and p are the proton, neutron and total densities, re- 
spectively. In the RPA, mRPA and ERPA calculations 
the one-body amplitudes x'^^, are defined using a large 
number of single-particle states including those in the 
continuum. We discretize the continuum states by con- 
fining the wavefunctions in a sphere with radius 15 fm 
and take all the single-particle states with ea < 50 MeV 
and ja < ll/2?i. As the residual interaction, we use the 
same simple force as that used in the ground-state cal- 
culation. Since the residual interaction is not consistent 
with the effective interaction used in the calculation of 
the single-particle states, it is necessary to reduce the 
strength of the residual interaction so that the spurious 
mode corresponding to the center-of-mass motion comes 
at zero excitation energy in RPA. We found the reduction 
factor / is 0.62. The p-h interaction in the second term 
on the right-hand side of Eq. Qj is multiplied with this /. 
To define the two-body amplitudes X^g^,^,, we use the 
small single-particle space consisting of the lsi/2, IP3/2, 
lPi/2, 1(^5/2, 2si/2, 14/2, 2p3/2, 2pi/2 and I/7/2 for both 
protons and neutrons. To reduce the number of the two- 
body amplitudes, we consider only the 2p-2h and 2h-2p 
components of Xj^„^,„, with \ea + e^ — e^' — e^'| < 60 
MeV. Since the single-particle space used for X^„^,„, is 
small, we use f — 1 for the matrix elements of the resid- 
ual interaction which couple to X^a^,0,. 

The occupation probabilities calculated in TDDM are 
shown in Table ID The deviation from the HF values 
{uaa^l or 0) is more than 10%, which means that the 
ground state of ^^O is a strongly correlated state. A 
recent shell-model calculation by Utsuno and Chiba [4I 
also gives a similar result for the ground state of ^^O. 
The correlation energy Ec in the ground state, which is 



TABLE I. Single-particle energies to 
bilities Uaa calculated in TDDM. 



and occupation proba- 





e. [MeV] 


^aa 


orbit 


proton neutron 


proton neutron 


lP3/2 
lPl/2 
1^5/2 
2Si/2 


-18.3 -21.9 

-12.3 -15.7 

-3.8 -7.1 

1.1 -1.6 


0.894 0.893 
0.868 0.865 
0.108 0.109 
0.019 0.021 
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FIG. 1. (Color online) Strength functions calculated in RPA 
(blue), SRPA (red), mRPA (grey) and ERPA (black) for the 
octupole excitation in ^®0. The distributions are smoothed 
with an artificial width F — 0.5 MeV. 



defined by E, = Ea/3a'/3'(a/3|«|a'/5')Ca'^'a/3/2, is -23.7 
MeV. A large portion of the correlation energy is com- 
pensated by the increase in the mean-field energy due 
to the fractional occupation of the single-particle states. 
The resulting energy gain due to the ground-state corre- 
lations, which is given by the total energy difference be- 
tween HF and TDDM, is 6.4 MeV, which is much smaller 
than \Er\ = 23.7 MeV. 
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FIG. 2. (Color online) Strength functions calculated in RPA 
(blue) and SRPA (red) for the isovector dipole excitation in 
^^0. The distributions are smoothed with an artificial width 
F = 0.5 MeV. 
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FIG. 3. Strength functions calculated in mRPA (grey) and 
ERPA (black) for the isovector dipole excitation in ^''O. The 
distributions are smoothed with an artificial width F = 0.5 
MeV. 



The strength functions for the isoscalar octupole ex- 
citation calculated in RPA (blue), SRPA (red), mRPA 
(grey) and ERPA (black) are shown in Fig. [T] The ex- 
citation operator used is r^Y^o^O). In mRPA and ERPA 
the fractional occupation of the 2si/2 states makes the 
transitions from these states to the continuum states 
possible, which gives some unphysical strength distribu- 
tions in very low-energy region. As will be shown below, 
the unphysical contributions could in principle be elim- 
inated by using the correlation matrix Ca'P'ap- How- 
ever, it is difficult to completely elimination the unphys- 
ical states. Therefore, we neglected the single-particle 
transitions from the 2si/2 states. In RPA the summed 
energy-weighted strength exhausts 97% of the sum rule. 
Other approximations give a similar value of the energy- 
weighted sum of the octupole strength. The excitation 
energy of the first 3" (3^;) state in RPA is 3.54 MeV, 
which is much smaller than the experimental data 6.13 
MeV [8|. In SRPA the state corresponding to the 3|f 
state becomes an imaginary solution due to the cou- 
pling to the 2p-2h configurations. The results in RPA 
and SRPA suggest that the aproaches based on the HE 
ground state cannot give a correct excitation energy of 
the 3j~ state. The excitation energy of the 3^ state in 
mRPA is 9.1 MeV, which is much larger than that in 
RPA. This is explained by the self-energy contributions 
(see Eq. (jj])) in the p-h pairs such as (lpi/2)~^ x 2^5/2 
and (1^3/2)^^ X 2(i5/2- The self-energy contributions in- 
crease the energy of the p-h pairs, reflecting the fact that 
the ground-state energy is lowered by the ground-state 
correlations. We found that the last four terms in Eq. 
@ , which describe the modifications of the p-h interac- 
tion, play a role in slightly increasing the attractive p-h 
correlations. In ERPA the 3j~ state is shifted downward 
to 6.56 MeV due to the coupling to the 2p-2h configu- 
rations. The B{E3) value of the 3" state at 6.56 MeV 
calculated in ERPA is 130 e^fm^, while the experimen- 
tal value is 204 ± 6 e^fm^ 1 81. The results in mRPA 



and ERPA demonstrate that both the ground-state cor- 
relations and the coupling to 2p-2h configurations play 
an important role in describing the properties of the 3^ 
state. 

We point out that the 3j~ state in RPA somewhat 
depends on the parameters of the Skyrme interactions. 
For example, it has been reported that the SkM* pa- 
rameter set gives the 3j" state at 6.06 MeV with 
B{E3) = 91.1 e'^m^ ^ and that a simple Skyrme 
type interaction with to = -1048 MeVfm^, is = 19150 
MeVfm^ and wq = 95 MeVfm^ gives the 37 state at 
6.05 MeV which exhausts 10.5% of EWSR [l|. We also 
found that the parameter set of the Skyrme III force with 
/ = 0.62, xo = and wq = 90 MeVfm^ gives the 3J" 
state at 6.12 MeV with B{E3) = 211 e^fm^ in RPA, 
which are close to the experimental data 6.13 MeV and 
204 ±6 e^fm^ Q. However, this force also induces strong 
ground state correlations and the 3j~ state calculated in 
ERPA comes at 8.26 MeV, which is much higher than 
the data. 

The strength functions for the isovector dipole excita- 
tion calculated in RPA (blue) and SRPA (red) are shown 
in Fig. [3 and those in mRPA (grey) and ERPA (black) in 
Fig. [3] The single-particle transitions from the partially 
occupied 2si/2 states are neglected in mRPA and ERPA. 
In RPA the summed energy-weighted strength exhausts 
84% of the dipole sum rule including the enhancement 
term, which is given by the ti and ^2 parameters of the 
Skyrme HI force and accounts 26% of the sum rule value. 
To increase the sum of the energy-weighted strength, we 
need to include the momentum dependent terms of the 
Skyme force in the residual interaction and also expand 
the single-particle space. The energy-weighted sums of 
the dipole strength in the other approximations give a 
value similar to that in RPA. Comparison of the result 
in SRPA with that in RPA shows that the main effect of 
the coupling to the 2p-2h configurations is to shift down- 
ward the RPA strength distribution. A similar downward 
shift of the dipole stren gth has been reported in large 
scale SRPA calculations [12, HJI • In mRPA the largest 
peak is upwardly shifted to 21.5 MeV from the position 
at 19.5 MeV in RPA. This upward shift is due to the 
self-energy contributions mainly in the (1^3/2)"^ — 1^5/2 
pairs. In the case of the isovector dipole excitation, the 
last four terms in Eq. ^ play a role in slightly reducing 
the repulsive p-h correlations. The increase in the dipole 
strength below 15 MeV in mRPA and ERPA, which is 
consistent with the experiment [lj|, is due to the partial 
occupation of the 1^5/2 states allowing such low-energy 
p-h transitions as 1^5/2 — > 2p3/2- The strength distribu- 
tion above 20 MeV in mRPA is about 2 MeV downwardly 
shifted in ERPA due to the coupling to the 2p-2h config- 
urations. However, there is little difference between the 
mRPA and ERPA distributions in the energy region be- 
low 17 MeV. This indicates that the dipole states in the 
low-energy region weakly couple to the 2p-2h configura- 
tions. In the case of the isovector dipole excitation the 
effects of the ground-state correlations play a role in in- 




FIG. 4. (Color online) Strength functions calculated in mRPA 
(black) and RPA (red) for the isovector dipole excitation in 
^^0. The partially occupied single-particle states obtained in 
TDDM are used in the RPA calculation. The correlation 
matrix Ca'p'ap is inserted only to the 2si/2 - '2.P3/2 and 2si/2 
- 2pi/2 pairs in the mRPA calculation. The distributions are 
smoothed with an artificial width F = 0.5 MeV. 



creasing the fragmentation of the dipole strength in low- 
energy region. However, in the giant dipole resonance 
(GDR) region it seems that the increase in the energy of 
the p-h pairs is compensated by a downward shift of the 
dipole strength due to the coupling to the 2p-2h configu- 
rations. The experimental photoabsorption cross section 
|14| shows a broader GDR distribution than the result in 
ERPA. More two-body configurations should be included 
to improve the ERPA result. 

As mentioned above, partial occupation of the particle 
states could cause some unphysical low-lying transitions 
possible. To show this, we performed an RPA calculation 
for the isovector dipole excitation using the partially oc- 
cupied single-particle states including the 2sii2 states. 
The result is shown in Fig. U with the red line. The 
strength seen below 5 MeV comes from the transitions 



from the 2si/2 states to the 2^3/2 and 2pi/2 states. The 
black line in Fig. S] shows the result in mRPA where 
the correlation matrix Ca' p'afi is included only in the 
2si/2 — 2p3/2 and 2si/2 — 2pi/2 p-h pairs. Since the norm 
matrix Si given by Eq. ([Sj for the 2si/2 - 2p3/2 and 
2si/2 - 2pi/2 pairs is quite small (w 0.02, see TableU]), the 
terms in Eq. ^ which contain Capa'P' can drastically 
shift the energies of these pairs. In fact the inclusion of 
the correlation matrix eliminates the strength in the low- 
energy region and slightly increases the strength distri- 
bution around 25 MeV, as seen in Fig. 2] When the cou- 
pling of the p-h pairs to 2p-2h configurations is included 
in ERPA, it becomes difficult to completely eliminate the 
unphysical low-energy transitions from the partial occu- 
pied 2si/2 states. This is the reason why we neglected 
such transitions in the niRPA and ERPA calculations 
shown above. 

In summary, the effects of the correlations in the 
ground state of "'^^O on the octupolc and dipole excita- 
tions were studied using the extended RPA. It was found 
that the ground-state correlation effects on the first 3~ 
state are significant: the first 3~ state is shifted upward 
due to the self-energy contributions. It was also found 
that the coupling to two particle - two hole states plays a 
role in producing the first 3^ at right excitation energy. 
In the case of the isovector dipole excitation, the effects 
of the ground-state correlations were found to increase 
the fragmentation of the dipole strength in low-energy 
region. However, the giant dipole resonance calculated 
in ERPA is little changed from that in RPA due to the 
competing effects: the increase in particle - hole energy 
and the coupling to two particle - two hole configurations. 
The mechanism to eliminate unphysical low-energy tran- 
sitions originating from partial occupation of the parti- 
cle states was also discussed. Our results demonstrate 
that the ground-state correlation effects in ^^O should 
be properly taken into account in the study of collective 
excitations. It is interesting to study these effects in other 
nuclei. 
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